Vortex localization in rotating clouds of bosons and fermions 
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Finite quantal systems at high angular momenta may exhibit vortex formation and localization. 
These phenomena occur independent of the statistics of the repulsively interacting particles, which 
may be of bosonic or fermionic nature. We analyze the relation between vortex localization and 
formation of stable Wigner molecules at high angular momenta in the view of particle-hole dual- 
ity. Trial wave functions for the vortex states and the corresponding fermion-boson relations are 
discussed. 
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PACS numbers: 03.75.Lm, 05.30.Fk, 73.21. La 
I. INTRODUCTION 

When small quantal systems are set rotating, vortices 
may form. These vortices are characterized by rotational 
flow of particle current around minima in the density 
distribution. In contrast to fluid mechanics, a vortex 
in a quantum system can be n-fold quantized, with the 
corresponding wavefunction having a (n-fold) zero at the 
vortex position and a phase changing by n2n on a path 
around this zero. 

The existence of a triangular vortex lattice in super- 
conductors was predicted by Abrikosov already in the 
50ies. In the mesoscopic regime, finite-size effects de- 
termine the symmetry of multivortex configurations. In 
small systems, a compromise must be found between the 
triangular lattice of the bulk, and the shape of the con- 
finement. It was predicted @,HHH 

that on a small su- 
perconducting disk, instead of the triangular Abrikosov 
lattice, the vortices may form geometric patterns that 
resemble those of classical point charges in a harmonic 
potential, i.e. small, finite-size Wigner crystals or so- 
called Wigner molecules 0, Q ■ This was experimentally 
confirmed very recently by Grigorieva et al. Q . 

However, the appearance of vortices and their localiza- 
tion is not limited to superconductivity or superfluidity. 
Another well-known example is a Bose-Einstein conden- 
sate (BEC) of (repulsive) alkali-metallic atoms (see, for 
example, Dalfovo et al. 9] and Legett or the book 
by Pethick and Smith [llj for reviews on BEC), which 
is set rotating. With increasing angular momentum, an 
ever larger number of vortices may penetrate the bosonic 
cloud of atoms [I1EG3. 

In BEC, much work focused on the Thomas- Fermi 
regime of strong coupling 0, 0, 0, 0, 0, 0] , where 
the kinetic energy dominates. In the dilute limit, how- 
ever, even though the interactions are weak they play the 
major role [U HJ H| . Within the Gross- 
Pitaevskii mean field approach, Butts and Rokshar [2l| 
and later Kavoulakis et al. [2J] found geometric vortex 
configurations between distinct values of angular momen- 
tum [271 l2Sj . In analogy to mesoscopic superconduct- 
ing disks, these vortex configurations again show selected 
symmetries corresponding to those of small Wigner crys- 



tallites |29j . 

Surprisingly, vortices may form independent of the 
statistics of the quantum particles. A striking identity 
between boson and fermion many-particle configurations 
leads to similar vortex states for bosons as well as for 
fermions, implying that vortex formation indeed is a uni- 
versal phenomenon in a rotating quantum system |30| ] . 

At extreme angular momenta, far beyond the point 
where the vortex lattice melts [3l| , bosonic 133, 133 . 134 
as well as the fermionic particles crystallize |34l 1.3-1 . l36| , 
approaching the classical limit j^l ■ 

The purpose of this article is to show how vortex for- 
mation and localization, as well as the crystallization of 
particles at the limit of extreme rotation for bosons and 
fermions are connected. Our analysis is largely based on 
the fact that vortex formation can most easily be ana- 
lyzed in terms of particle-hole duality, which was found 
to hold for the bosonic as well as the fermionic case [2£| • 

This paper is organized as follows: In section II, af- 
ter describing the model and methods, we briefly discuss 
the particle-hole dualism for the fermion vortex states. 
Results of the direct numerical diagonalization of the 
many-body Hamiltonian follow in III, where we discuss 
the localization for bosons and fermions at high angular 
momenta in terms of the regular oscillations in the yrast 
line. We apply similar arguments to explain the localiza- 
tion of vortices, both in the bosonic and fermionic limit. 
For fermions, this fact is further illustrated making use of 
particle-hole duality. Trial wave functions for the vortex 
states and the corresponding fermion-boson relations are 
discussed in section IV. 



II. THE MODEL 

A. The Hamiltonian 

Let us now consider interacting particles confined by 
a two-dimensional harmonic trap. Ignoring the spin de- 
gree of freedom, these particles can either be spinless 
bosons (as for example, bosonic atoms in the same hy- 
perfine state), or polarized spin-1/2 fermions (say, elec- 
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trons). The many-particle Hamiltonian is simply 
n 2 n N N 

i i i<j 

where N is the number of particles with mass fi, lu the 
oscillation frequency of the confining potential, and v(r) 
a repulsive two-body interaction. For a (dilute) gas of 
spinlcss, bosonic atoms, often a contact interaction of 
the form v(r) = ^-Uq8{t) is used, where t/ then depends 
on the scattering length for the atom-atom collisions [24| . 
For spinless fermions, however, this interaction only re- 
covers the non-interacting case due to the Pauli princi- 
ple. To ease the direct comparison between the boson 
and fermion spectra, in both cases we thus consider the 
usual long-range Coulomb form, v(r) — e 2 /4ire r. For 
angular momenta well below the fractional quantum Hall 
regime [2(J |2^, |38| , the boson spectra, calculated with ei- 
ther the short-range or long-range Coulomb interaction, 
in fact show a remarkable similarity, as earlier demon- 
strated by Toreblad et a l. 1301 . We refer to the work 
by Hussein et al. [U ES Mfand Vorov et al. ;27] for a 
more general discussion of the universality of repulsive 
interactions in the bosonic case. 

In what follows, we are mainly interested in the gen- 
eral structure of the eigenstates of Eq. 2] at high angular 
momentum M. We follow the nuclear physics tradition 
to study the so-called yrast line, i.e. the lowest-energy 
states as a function of M, and the correspon ding low- 
lying excitations (the so-called yrast spectrum |42l l43|l. 
Similarly, magnetic fields are not included explicitly: ne- 
glecting the Zeeman term, the only effect a magnetic field 
has, is to set the system rotating. 

B. Restriction to the Lowest Landau Level 

In the absence of interactions, the single-particle en- 
ergies of the two-dimensional harmonic oscillator are 
e = huj(2n + \m\ + 1), where n is the radial quantum 
number, and m the single-particle angular momentum. 
The energy level structure is schematically sketched in 
Figure At very large total angular momentum of the 
non-interacting many-particle system, the lowest-energy 
state is characterized by quantum numbers n — 0, and m 
being zero or having the same sign as the total angular 
momentum, M. Measured relative to the non-interacting 
ground state, the total energy of the lowest state at a 
given M thus equals hu>(M+l). This single-particle basis 
is identical to the so-called Lowest Landau Level (LLL) 
at strong magnetic fields. In this subspace, a configu- 
ration can be denoted by the Fock state |rio^i^2 ' ' ' n s), 
where n, is the occupation number for the single-particle 
state with angular momentum i, and s the largest single- 
particle angular momentum included in the basis. For 
spinless fermions, ni is or 1 and for bosons it is an in- 
teger. As the angular momentum M is a good quantum 
number, we have the restriction irii — M . 
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FIG. 1: Single-particle energy levels of a two-dimensional har- 
monic oscillator. The black bullets show the Maximum Den- 
sity Droplet (MMD) configuration for six particles, and the 
open circles the M = ground state for noninteracting, polar- 
ized fermions. The lowest-energy levels for each m form the 
Lowest Landau Level (LLL), as indicated by the bold lines. 



For fermions, the smallest angular momentum that can 
be built in the LLL is that of the so called Maximum Den- 
sity Droplet 0], Mmdd = N(N - l)/2, where due to 
the Pauli principle the lowest-energy states with succes- 
sive single-particle angular momenta are occupied, gen- 
crating a compact configuration 

| 111... Ill 000...) (2) 

N 

(see Fig. Pi . 

To see the similarity with quantum Hall systems [45| . 
one may define a filling factor v by the ratio between 
Mmdd and the actual angular momentum M, 

v = N(N — 1)/(2M) . (3) 

Clearly, the MDD corresponds to the Laughlin quan- 
tum Hall state at filling factor one in the bulk ^(|. For 
bosons, the state with smallest possible angular momen- 
tum within the LLL is the non-rotating ground state, 
where all particles occupy the single quantum state with 
n = and m = 0, i.e., | A^OOO . . . }. 

C. Numerical diagonalization of the Hamiltonian 

Including now the interactions between the particles, 
we solve the full many-particle Hamiltonian exactly by 
numerical diagonalization, in order to obtain the ground 
state and the low-lying excited states at a given angu- 
lar momentum M. In this so-called configuration in- 
teraction (CI) method, the many-particle state is de- 
scribed as a linear combination of Slater determinants, 
I*) = J2Ca\£-a), which for given statistics are defined 
by the combinatorics of the single-particle basis. Here, it 
is natural to use the basis of the Harmonic oscillator (i.e. 
that of the external trap) . If the strength of the interpar- 
ticle interaction is much smaller than the single-particle 
excitation energy hu, at large M the most important con- 
figurations indeed consist of single-particle states of the 
LLL. In the limit e 2 / 1 Aitet)r s hujQ — > (or v /huJo — > 0), 
where r s is the electron density parameter, uq = l/Vr 2 , 
configurations with particles at higher Landau levels have 
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diminished weight C a , and the LLL approaches the exact 
result. (For a discussion that includes Landau level mix- 
ing in the small- N limit, see the recent work by Stopa et 
al. H3 , as well as Giiglii et al. 0> @| > considering quan- 
tum dots at high magnetic fields. A similar study was 
carried out for trapped bosons with contact interactions 
by Morris and Feder [50|). 

Increasing the angular momentum beyond the MDD in 
the fermionic case, M > Mmdd, or beyond the "conden- 
sate" , M > 0, in the bosonic case, there are many ways to 
distribute M quanta of angular momentum on the given 
number of particles, N. For large values of N and M, 
we thus have to restrict the basis by setting a smaller 
value s for the largest single-particle angular momen- 
tum. However, the most dominant configurations appear 
rather compact, i.e. most of the occupied single-particle 
states are close in angular momentum m. Thus, accurate 
results can still be obtained in a restricted Hilbert space. 

In the LLL the single-particle basis states in polar co- 
ordinates (r, 4>) are 

4> m (r,<t>)=A m r m e-^' i e im +, (4) 

where A m is a normalization factor. 

The many-particle Hamiltonian, Eq. ft}, can be writ- 
ten, apart from a constant, as 

H = '^2m l hujcjc i + ^ VijkicfcfckCi, (5) 

i i,j,k,l 

where 

V m = J I dvdr'v(\r-r'\)^(r)tp*(r')Mr)Mr'), (6) 

and v is the interparticle interaction (see above). Since 
the total angular momentum is fixed, M = X) m *> the 
diagonal part gives an energy TiujM for all configurations. 
In effect, we thus only diagonalize the interaction part of 
the Hamiltonian. This simplification is caused by the 
restriction of the basis to the lowest Landau level and 
holds for bosons as well as for fermions. It is natural to 
present the energies of the interaction part only (which 
trivially scales with its strength) , as done for most of the 
results shown below. 

The interaction matrix elements are computed numer- 
ically. In the case of Coulomb interaction we used the 
technique suggested by Stone et al. 51]. For numerical 
diagonalization of the Hamiltonian matrix, we applied 
the Lanczos method l52l. 



D. Particle- hole duality 

A Fock state can be described either by particles, or 
equivalently, by holes. For spinless fermions, in prac- 
tice one simply can replace the zeros with ones, and vice 
versa the ones with zeros. For example, the configuration 
|1110011111) is then changed to 0001100000). One cor- 
respondingly defines the creation (annihilation) operator 



of a hole as d + = c (d = c + ) and writes the Hamiltonian 
Eq. JSJ in terms of these: 

H = «^o(l ~ d+di) + 2^2 (Vijij ~ V l]Jl )d+d l 

i i,j 

+ 22 Vijkidfdfrdjdi + constant. (7) 

i,j,k,l 

Naturally, the solution of this Hamiltonian leads to an 
equivalent result, and requires the same numerical effort. 

The interactions between the holes here are the same 
as the interactions between the particles. However, it 
is important to note that the second term is nonzero. 
This means that the holes effectively do not move in a 
harmonic confinement. This has two important conse- 
quences: the holes (i) do not have a pure center of mass 
excitation, and (ii) not necessarily the same localization 
geometry as the particles. 

This particle-hole dualism, as earlier described in 
Ref. j2^, allows us to identify the holes as vortices in 
the Fermi sea and to evaluate the correlation functions 
between these vortices. Furthermore, in the case of a few 
holes in the Fermi sea, the energy spectrum is dominated 
by these holes and can be understood by diagonalizing 
only the interaction part of the hole Hamiltonian (i.e. the 
third term in Eq. Q). 

III. RESULTS 

A. Localization of particles at extreme angular 
momenta 

Let us first analyze the many-body quantum spectra 
for small particle numbers. (This will later turn useful 
when studying the spectral properties of many-particle 
systems in terms of localized vortices) . 

Figure [21 compares the spectra for three bosons (red 
crosses, +) and three fermions (blue crosses, x) at ex- 
treme angular momenta, M > ZMmdd- (Only the non- 
trivial interaction part of the total energy is shown, here 
for the ten lowest states at fixed M). Naturally, the in- 
teraction energy goes down with increasing angular mo- 
mentum for repulsive interactions. Note that for each 
new energy at a given value M' , the same energy is re- 
produced at all angular momenta M > M' . These are 
simply the center of mass excitations characteristic for 
harmonic confinement |5^]. With increasing M, the bo- 
son and fermion spectra become remarkably similar. This 
does not only hold for the yrast line, but also for the low- 
lying excitations. 

The yrast line (drawn in Fig.[3as a red line to guide the 
eye) shows clear oscillations with a period of AM = 3. 
The inset shows the corresponding spectrum for N = 2, 
with a similar oscillation in the yrast states, alternating 
between even and odd angular momenta, respectively. 
Further increasing the number of particles, we show in 
Fig. |3| the yrast spectrum for N — 6 bosons. Again, the 
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FIG. 2: Interaction energy of a three-particle system, N = 3, 
for bosons (red crosses, +) and fermions (blue crosses, x), as 
a function of the angular momentum, for angular momenta 
AI > 3Mmdd ■ (The inset to the right shows the same for 
N = 2). The red (blue) line is the yrast line. The insets 
to the left show the boson and fermion correlation functions, 
here taken at the cusp states with M = 21 as an example. 
The white dot indicates the reference point. The particles 
appear localized at the corners of a triangle, with one particle 
placed at the reference point. 



boson and fermion spectra (not plotted here) are very 
similar in the large- M limit. 

Similar to the results shown for N — 3 in Fig. above, 
the boson spectrum in Fig. [3] shows pronounced cusps 
with a characteristic gap to the low- lying excitations. At 
smaller M, oscillation periods with AM = 5 and AM = 
6 compete, with cusps at M = 18,20,24,25,30,36,40.... 
Beyond M = 50 the cusps appear regularly with AM = 
6. This is very clearly seen in the inset to Fig. [21 where 
a third order polynomial in M was subtracted from the 
energies, changing the slope of the yrast line in order to 
make its oscillations more visible. 

These cusp states and their periodicities in M are a 
phenomenon well known to occur in quantum dots at 
very strong magnetic fields, which have been investigated 
extensively in the literature [45|]. Here, the regular oscil- 
lations of ground state energy as a function of magnetic 
field were found to originate from r igid rotations of the 
classical electron configurations [5J, [55| ■ While for par- 
ticle numbers up to N — 5, the Wigner crystallites form 
simple polygons (from triangle to pentagon), the case 
N = 6 is more complicated due to the interplay between 
two stable classical geometries, a pentagon with one par- 
ticle at the center, (1, 5), and a hexagon, (0, 6) |32l |. 

These simple geometries in fact easily explain the pe- 
riodic oscillations in the above many-body spectra. The 
rigid rotation of the 'molecule' of localized particles is 
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FIG. 3: Interaction energy of N = 6 bosons as a function of 
the angular momentum. The inset shows the yrast line with 
a smooth function of angular momentum (3rd order polyno- 
mial) subtracted from the energies, in order to make the oscil- 
lations more visible. The large-M limit is dominated by a reg- 
ular oscillation with AM = 5. The pair correlation functions 
to the left clearly demonstrate localization in Wigner molecule 
geometries at high angular momenta. While at smaller M- 
values, the (1,5) and the (0,6) configurations compete, at 
extreme angular momenta fivefold symmetry dominates. 

only possible at angular momenta allowed by the under- 
lying symmetry group |34j . For example, if the electrons 
localize for example in an equilateral triangle, the three- 
fold symmetry leads to period of three. The five-fold 
symmetry of the (1,5) ground state of the six-electron 
molecule yields a period of five EE EE E3- Cor- 
respondingly, the period AM — 6 is explained by the 
hexagonal structure, (0,6). 

A related periodicity of the yrast energy as a function 
of the angular momentum is obtained also for other par- 
ticle numbers, and holds even if higher Landau levels are 
included in the basis set 03 EE Ell ■ (Similar agru- 
ments were successfully applied earlier also to analyze 
the rotational spectrum of quasi-one-dimensional quan- 
tum rings, see EEE3)- 

The particle localization in the internal coordinates of 
the quantum system can be demonstrated further by us- 
ing a rotating frame |55| . or by studying the pair corre- 
lations defined as 

ff (r,r') = <*|n(r)n(r')|*> , (8) 

where n is the density operator. The pair correlation 
function g describes the probability to find a particle at 
r' when another particle is in the reference point r. In 
a finite system, the pair correlation function is a func- 
tion of two variables, and often also called "conditional 
probability". In Fig. |3 (le.fi) we show the pair correla- 
tions for the three-particle system, both in the bosonic 
and fermionic case, at M = 21. In both cases, the parti- 
cles appear localized in a triangular geometry, with two 
particles showing clear maxima opposite to the reference 
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point (indicated by the white dot). A very similar result 
is obtained for the six-particle system, where at smaller 
M-values, the (1,5) and the (0,6) configurations com- 
pete. At extreme angular momenta, however, as here 
shown for M = 60, the fivefold symmetry (1,5) domi- 
nates. These pair correlations clearly demonstrate that 
at high angular momenta, bosons as well as fermions lo- 
calize in Wigner molecule geometries, in excellent agree- 
ment with the earlier results by Manninen et al. [3^ | , and 
their subsequent confirmation by Romanovsky et al. |33| . 



B. Localization of Vortices 

For a harmonically trapped cloud of bosons that is set 
moderately rotating, geometric vortex configurations are 
known to appear between distinct values of angular mo- 
mentum. A single vortex at the center of the system is 
formed when M = N, two vortices at about M = 1.7N, 
three vortices when M — 2. IN and N vortices before 
M = 3N. It was shown earlier that these mean-field re- 
sults within the Gross-Pitaevskii scheme [^J, O0| emerge 
as the correct leading order approximation to exact cal- 
culations in the same subspace [59^. 

The many-body spectra for 20 and 40 bosons inter- 
acting by the Coulomb force, are compared in Figure 01 
displaying the yrast line and up to ten low-lying excita- 
tions. The horizontal axis is now the angular momen- 
tum per particle, / = M/N. As previously, a second- 
order polynomial in M was subtracted to emphasize the 
salient oscillations of the yrast line. For these large par- 
ticle numbers, the Hilbcrt space was truncated such that 
only single-particle states with m < 10 were included, 
still giving well-converged results in the bosonic case. 

Both spectra show pronounced oscillations with con- 
secutive periods of AM = 2,3 and 4, expanding through 
the above mentioned intervals of subsequent vortex entry. 
While a single vortex is formed at I = 1, oscillations with 
AM = 2 coincide with the two- vortex solution appearing 
between 1.7 < I < 2.1, three vortices beyond I = 2.1, etc. 
These results are in excellent agreement with the results 
of the Gross-Pitaevskii approximation [24| . 

The regular oscillations in fact appear very similar to 
the cusp states at extreme angular momentum, as dis- 
cussed for small particle numbers in Figs. Gland EUabove. 
However, in contrast to localization in Wigner molecules, 
these oscillation are now caused by two, three, or four 
localized vortices, respectively. (The 2- and 3-vortex ge- 
ometries are schematically sketched in Figs. 0] and Elbe- 
low). 

We notice that the spectra look qualitatively very sim- 
ilar to that for 20 fermions shown in Fig. There is, 
however, a marked quantitative difference in the angu- 
lar momentum where the vortices appear. In the case 
of bosons, the value of angular momentum per particle 
M/N where the second, third, etc. vortex enters the 
cloud, seems to be independent of the number of bosons 
in the system. (Naturally the number of particles has to 
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FIG. 4: Many-particle energy levels for 20 (lower panel) and 
40 bosons (upper panel) as a function of angular momentum 
per particle, I = L/M. In both cases, a second-order poly- 
nomial in M was subtracted from the energies to make the 
yrast line nearly horizontal. 



be much larger than the number of vortices) . If the same 
systematics would hold for fermions, we should expect 
the same small number of vortices to appear when the 
angular momentum is increased beyond the MDD. This, 
however, does not hold for N > 12. The vortex systemat- 
ics then changes [6(3, 0, as compared to the bosonic 
case. The increase of angular momentum per particle, 
(M — Mmdd)/N, for the appearance of the second vor- 
tex depends on N, as illustrated in Table |U 

The origin of the periodic oscillations can be under- 
stood by studying the corresponding many-particle prob- 
lem of holes, as briefly discussed in section llTDl above. 

For this purpose, let us now study in more detail the 
fermion case. The excitations from the MDD to higher 
angular momenta create holes in the filled Fermi sea. 
These holes are associated with the existence of vor- 
tices in the system 29] . In the single- vortex region, the 
most important configuration has a single hole at angu- 
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TABLE I: Angular momenta where the vortices appear for 
boson and fermion systems of different sizes. N is the number 
of particles, M the angular momentum and M = M—Mmdd 



vortex 


[M/N] for bosons 
N = 8 N = 20 N 


= 40 


[M /N] for fermions 
N = 8 N = 20 


2nd 


1.7 


1.6 


1.6 


1.7 


0.9 


3rd 


2.2 


2.1 


2.1 


2.2 


1.4 


4th 




2.8 


2.8 
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FIG. 5: Many-particle energy levels for 20 fermions as a func- 
tion of angular momentum per particle relative to the MDD, 
I = (M — Mmdd)/N. As in Fig. 2] above, a second-order 
polyonomial was subtracted from the yrast line. 



lar momentum Mh = Mmdd + N — M. When M is 
increased from the MDD the single vortex enters the sys- 
tem from the surface and reaches the center at angular 
momentum M = Mmdd + N. Then, Mh = 0, i.e. the 
Fock state with the largest amplitude has the structure 
| 0111 . . . 1111111000 . . .}. This state appears in the ex- 
citation spectrum and terminates the band [(5lJ : Within 
the lowest Landau level, no further quanta of angular mo- 
mentum can be added, unless a second hole penetrates 
the electron droplet. Note that the single hole has a sim- 
ple single-particle wave function, as in Eq. i.e. the 
hole (or vortex) is either localized at the center, or delo- 
calized on a ring. 

The single vortex does not reach the origin before it is 
energetically more favorable to create two vortices that 
are closer to the surface. Considering the state for holes, 
the angular momentum is written as 



M h 



(N + h)(N + h-l) 



M, 



(9) 



where h is the number of holes and M the angular mo- 
mentum of the fermion system. For example, the 20- 
particle state with angular momentum M = Mmdd + 
24 = 214 corresponds to a two-hole state with Mh = 17. 
The two-hole state at such high angular momentum is 
strongly correlated, with holes localizing in a narrow ring, 
as shown in Fig. below. 

Similar to the particle-particle (pp) correlations, we 



can define a hole-hole (hh) correlation function 
gh h (T,T') = (V\h(r)h(r')\*), 



(10) 



where h(r) = Y^ij' , Pi( r ) t l ; j( r )^t^'j 1S ^ ne n0 ^ e density 
operator. 

The pair correlation between the holes clearly shows 
that in the internal frame, the holes indeed are local- 
ized. Fig. [|)| shows the that the density profile of the two 
holes forms a narrow ring. At the same radius, the hh- 
correlation shows a very well localized maximum, corre- 
sponding to a hole (vortex) localized at the opposite side 
of the "reference hole" . 

N=20 
M- 214 




density profile 
2 holes 



hole-hole correlation 



FIG. 6: Density profile of two holes at angular momentum 
Mh = 17 (a) and the corresponding hole-hole correlations 
(b). The hole density and pair correlation was determined 
from the calculation for 20 electrons with angular momentum 
M — 214 and the minimum single-particle basis allowing that 
angular momentum. 

Considering the two-hole system as a rigid "molecule" 
of two vortices, the oscillations in the energy spectrum 
consequently are understood as resulting from the rota- 
tional states: Like for a two-atomic molecule, only every 
second angular momentum is allowed in the rotational 
spectrum |63| ■ At angular momenta in between, the rota- 
tional state must be occupied by a vibrational excitation. 
Thus, the state has higher energy. Similarly, this holds 
for the periodic oscillations of the two-particle yrast line 
for bosons as well as fermions, as displayed in the inset 
to Fig. [3 above. 

Increasing the angular momentum, the whole system 
expands, while the vortices move closer to the center. 
Beyond a certain angular momentum, however, the re- 
pulsion between the vortices becomes so strong that it is 
energetically favorable to form a larger ring with three 
vortices instead. In this region, the spectrum (Fig. |3J) 
then shows an oscillation period with AM = 3, caused 
by the rigid rotation of the triangle of three vortices. Fur- 
ther increasing the angular momentum adds additional 
vortices, with the oscillation period increased by one each 
time a new vortex enters the system. 

Let us now study in more detail the three-vortex re- 
gion, and compare the spectrum calculated for only three 



7 



fermions, as shown in Fig- El with that for 20 fermions, as 
displayed in the lower panel of Fig. [7] (For a quantitative 
comparison of the spectra, we first subtracted from the 
20-fermion spectrum a linear function 0.9M and then 
plot it as a function of the corresponding hole angular 
momentum as defined in Eq. ijffil). 



Spectrum for 
20 electrons 



20 22 24 26 23 30 32 

Angular Momentum M | Ti ] 



Spectrum for 
3 holes 



24 26 23 

Angular Momentum M [ n ] 



FIG. 7: (a) Many-particle spectrum for three holes, as de- 
termined from the 20-fermion spectrum (see text). A linear 
function 0.9M was subtracted. The hole spectrum compares 
well with that of three particles, shown in (b) for Coulomb 
interactions. The horizontal lines in (b) connect states which 
are center-of-mass excitations. To guide the eye, correspond- 
ing lines are plotted in (a). 

The spectra are strikingly similar, indicating that in- 
deed the spectrum of 20 particles (in this region) is de- 
termined by the three localized holes. The lines in Fig.0 
indicate the center of mass excitations for three parti- 
cles. In the 20-particle spectrum, these states are not 
pure center of mass excitations of holes, since the true 
Hamiltonian for holes, Eq. (JJJ, has additional terms with 
anharmonic corrections to the confining potential. 

The observation that these quantum states show much 
similarity in their internal structure is further clarified 
by comparing their pair correlations. Fig. [S] shows the pp 
correlations at three points along the yrast line. When 
the yrast line has a kink, as here for example at Mh = 15 
or M = 239, the correlation function clearly shows local- 
ization of holes in a triangular geometry. At Mh — 16 the 



three-particle system originated from the Mh = 15 state 
by a center of mass excitation, and at Mh — 17 it has 
an internal vibrational mode yj, |23 • The effect of these 
excitations is clearly seen in the correlation functions. 
The correspondence between hh and pp correlations for 
N = 3 and N = 20, is not perfect, however, since in 
the 20-particle system the center of mass excitation (for 
holes) is mixed with the vibrational excitation. 



FIG. 8: The lower panel shows contours of the /i/i-correlations 
for three-vortex states of a 20-fermion system (Coulomb) at 
angular momenta M = 230, M = 231 and M = 232. The 
upper panel shows the corresponding pp-correlations for three 
particles, at M = 15, 16 and 17. 

The three-particle and the three-hole states for N = 20 
can be compared directly by calculating their overlaps. 
In doing so, the electron calculation has to be restricted 
to the minimum Hilbert space having only three holes. 
As an example, we calculated the overlap matrix for the 
four lowest energy states at angular momentum M = 24 
for the holes (the corresponding 20 electron angular mo- 
mentum is 229). The result is shown in Fig. 03 The over- 
lap between the ground states of the two calculations is 
large, showing again that their internal structure is simi- 
lar. The first and second excited states seem to be mixed, 
while the third excited state is again quite similar. This 
result is in agreement with the pair correlation functions 
(Fig. |SJ and the spectra (Fig. 0) which show that a pure 
center of mass excitation of holes (first excited state here) 
does not exist in the 20-particle spectrum. 

For small particle numbers, N ~ 10, the pair corre- 
lations cannot display the vortex structure properly. In 
this case, weak signals of vortices can be seen only if 
the reference point is clearly outside the density distri- 
bution j^l, since then the exchange-correlation hole is 
distributed rather evenly over the remaining dot area. 

The situation becomes different, when the particle 
number is significantly increased. Figure 1101 shows the 
electron-electron pair correlation for 36 electrons. In the 
MDD, only the exchange hole at the reference point is 
clearly visible. In the case of M = 706, in addition 
to the exchange hole we see four wide minima, which 
we interprete as four localized vortices. The reference 
point is chosen at a radius where the density has a mini- 
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FIG. 9: Overlaps between 3-hole states from two different cal- 
culations for M h = 24 (M = 22). GS, ESI, ES2 and ES3 refer 
to the ground state, and excited states of the three particle 
calculation. The columns show the overlaps with the ground 
state (blue), and the first three excited states (red, yellow 
and green) of the hole state derived from the 20 electron- 
calculation. 



mum, i.e. at the expected radius of the vortex ring. The 
/i/i-correlation, also shown in Fig. ^3 shows three pro- 
nounced maxima, which are consistent with the minima 
in the pp-correlation. This comparison demonstrates, 
that the vortex localization is more clearly seen in the 
/i/i-correlations, than in the ee-correlations. The lowest 
panel in Fig. llOl shows the corresponding pair correlations 
in the bosonic case. Also in this case, we see four minima, 
corresponding to the four vortices, but naturally, not the 
exchange hole. 

Regarding the analysis of the corresponding Fock 
states in the bosonic case, it is not as straightforward 
to relate the vortices to holes in the similar fashion as 
for spinless fermions. Nevertheless, we can transform the 
boson wave function to a fermion wave function using the 
method explained in Ref. |30| and then plot the vortex- 
vortex correlation for the corresponding fermion system. 
Figure ^2 shows a pair correlation function constructed 
in this way, with the localization of the second vortex 
opposite to the reference point. 



Comparison between boson and fermion wave 
functions 



Comparing the energy spectra and the pair correlation 
functions, the vortex formation for bosons and fermions 
appears surprisingly similar. It is now interesting to see 
how far this similarity is reflected in the detailed struc- 
ture of the many-particle states. For this purpose, we 
need to compare the fermion wave function with angular 
momentum M — Mm unto the boson wave function with 
M. To this end, we should multiply the boson wave func- 
tion with the determinant Y[( z i ~ z j)- Here, we instead 



N=36 



Fermions 

particle-particle 
correlations 



hole-hole 
correlations 



MDD ™ M=70^^^ 



Bosons 




particle-particle 




correlations 






M=104 ^ 



FIG. 10: Pair correlation functions calculated for 36 electrons. 
The upper panel shows the electron-electron correlations for 
the MMD (left), for particles at M = 706 showing four vor- 
tices (right), and for holes at the same angular momentum 
(lower right). The lower panel shows the corresponding cor- 
relation function for a bosonic four-vortex state at angular 
momentum M = 104. (Note the absence of the exchange hole 
in the bosonic case.) 



Vortex-vortex correlation 




Bosons 



Fermions 



FIG. 11: Vortex- vortex correlation function of the two- vortex 
state for 20 bosons with M = 34, M/N = 1.7, left, and for 
20 fermions with M = Mmdd + 34 = 224. The fermion state 
is the second excited state for that angular momentum. For 
bosons the correlation function is determined by transforming 
the boson state to a fermions state and plotting the hole-hole 
correlation. 



use a simpler mapping based on the one-to-one corre- 
spondence between the boson and fermion configurations, 
as described earlier by Toreblad et al. |2jJ. The bose 
"condensate", |iV0000- • ■), corresponds to the MDD in 
the fermion case. Other configurations can be obtained 
as single-particle excitations out of the condensate. Ta- 
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ble ITT1 shows a few examples of these configurations, here 
for the simple case of six particles. Figure H"2l compares 



TABLE II: Examples of corresponding configurations for six 
fermions and bosons. M is the angular momentum of the 
configuration. For six fermions the angular momentum of the 
MDD is 15. 



A-^fcrmion 


fermion state 


■^boson 


boson state 


15+0 


11111110000) 





1600000} 


15+4 


11101111000} 


4 


1240000} 


15+4 


11111001100} 


4 


1402000} 


15+15 


11000111110} 





100500} 



the amplitudes of the most important configurations for 
the two-vortex state of 20 bosons with M — 34 to those 
of 20 fermions with M — 224. For bosons, the state 
is at the yrast line, while for fermions it is the second 
excited state. The yellow columns show the cumulative 
overlap between the boson and fermion states. The con- 
figurations are shown at the right, using the notation for 
bosons. The figure shows that for both kinds of particles, 
the same configurations are important. The actual am- 
plitudes differ, but the qualitative similarity of the states 
is clearly seen, especially regarding the signs of the differ- 
ent terms contributing to the linear combination of Slater 
determinants in the many-body state. 
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symmetry. Thus, only the radial density distribution is 
relevant. In the LLL the density can be simply deter- 
mined from the occupancies of the single-particle levels 



n(r) 



E 



Cal 2 E" 

\i> m {r,0)\ 



(11) 



where C a is the amplitude of the configuration a, n am 
the occupation number of the single particle state m 
in configuration cv, and ?p m the single-particle state of 
Eq. We note that the same information is contained 
in the occupancy of the single-particle angular momen- 
tum states: 



D(m) = ^ \C a 



(12) 



The restriction to the LLL makes it simple to determine 
the current density 



je{r) = 53 \ c u\ 2 ^2n mn — \ip m (r 7 9)\ 



m 
r 



(13) 



and the velocity field jg(r)/n(r). Naturally, the current 
density as well as the velocity field have the circular sym- 
metry of the underlying Hamiltonian. 

We can also use the vorticity of the velocity field 



V x 



Mr) 
n(r) 



ee = 



0_ 

Or 



Mr) 

n(r) 



(14) 



to give information about the vortex formation. If the 
vortex is at the center of the trap, naturally the density 
of the single- vortex state is rotationally invariant. How- 
ever, if one or more vortices are off-center, their charac- 
terization is not as straight-forward. 

For example, the two-vortex solution, as displayed in 
Fig. 1131 shows a clear maximum in vorticity where the 
radial density profile has a local minimum. The pro- 
files of the density and vorticity for bosons and fermions, 
respectively, appear rather similar, with the main differ- 
ence that the fermionic cloud extends to larger radius as 
a simple consequence of the Pauli principle. 



FIG. 12: Comparison of the coefficients of the most important 
configurations of the two- vortex state of 20 bosons with M = 
34 to those of 20 fermions with M = 224. In the case of 
fermions the state is the second excited state for the chosen 
angular momentum. For fermions, the amplitudes are shown 
in blue, and for bosons in red. The yellow columns show 
the cumulative overlap between the boson and fermion states. 
The configurations are shown at the right, using the notation 
for bosons. 



D. Density and vorticity 

Any exact solution of interacting particles in circularly 
symmetric potential must have a density with circular 



IV. TRIAL WAVE FUNCTIONS FOR VORTEX 
STATES AND FERMION-BOSON RELATIONS 

The fermion wave function for the maximum density 
droplet (MDD) is a single determinant, see Eq. where 
the number of singly occupied states equals the num- 
ber of particles. Using the single-particle states of the 
LLL, Eq. now in terms of the complex coordinates 
Zj = Xj + h/j , this takes the form of the Laughlin wave 
function |46j which, apart from the normalization, is writ- 
ten as 

= n^-^) 9 exp(-53|z fc | 2 /4j, (15) 



10 



FERMIONS 

N=20, L=224 

2nd excited state 2 vortices 



>. 08 

1 

Kg o.a 
Q 

o 

-0.3 





\ density 


vorticity 





2 4 8 8 

Radius r [a.u.] 

BOSONS 
N=20 L=34 

ground state, 2 vortices 



10 



c 
G 




2 4 S 

Radius r [a.u.] 

FIG. 13: Density (red lines) and vorticity {green lines) in 
a fermionic (upper panel) and bosonic (lower panel) system, 
respectively. For N = 20 particles, the fermionic two-vortex 
state at M = 190 + 34 is compared to its bosonic counterpart 
at M = 34. 



where q = 1 for MDD and in general an odd inte- 
ger for fermions. For the MDD, we use the notation 
<Ff = When the angular momentum is in- 

creased, there is no simple analytic formula for the exact 
state. However, it has been shown that at SMmdd, the 
next Laughlin state with q = 3 corresponding to the fill- 
ing factor v = 1/3 for the fractional quantum Hall effect 
is a good approximation to the exact wave function |32| . 
For angular momenta in between the Laughlin states, no 
simple analytic expressions for the wave function exist. 

Rather than providing accurate analytic estimates for 
the many-body wave function, our aim here only is to 
reach a physical understanding of vortex formation and 
crystallization at high angular momenta. To this end, 



let us first consider the simplest case of a single vortex. 
Bertsch and Papenbrock 23] have shown that for weakly 
interacting bosons, a single vortex at the center of the 
condensate can be described as 



N 



*f« = Y[( z i - z o) exp 




(16) 



where zq = {z\ + Z2 + ■ • ■ zjf)/N is the center-of-mass 
coordinate. We can use the same Ansatz for fermions 
and write 133 



JY 



F 

MDD- 



(17) 



Note that the only difference between the boson and 
fermion states is the additional product IIi<j( z » — z j) 
which makes the fermion state antisymmetric. We have 
tested this approximation in the case of six electrons. 
In this case, the overlap between the trial wave func- 
tion Eq. (fT7|) and the exact result (restricted to LLL) for 
Coulomb interaction is 98.5 %. Even when the restric- 
tion to LLL is abandoned and higher Landau levels are 
included, the overlap was found [33 to be 90% . 

The most important configuration of the Bertsch- 
Papenbrock state Eq. i|16|) f° r > sa y Slx > fermions is 

01111110000), i.e. one electron is removed from the 
angular momentum m — state and lifted to the first 
empty angular momentum state. We notice that this is 
the smallest angular momentum state that can have an 
empty state at m = 0. For six electrons, Mmdd + N — 
15 + 6 = 21. For smaller values of M the empty 
state is at a higher single-particle state. For example, 
for M = Mmdd + 1 the only possible configuration is 

11111010000) and the corresponding wave function is 
x ^mdd- When the angular momentum is increased, 
the hole in the Fermi sea moves to lower angular momen- 
tum until it reaches the origin. We will call these states 
as one-vortex states. When the hole is not at m = the 
vortex is not localized at the origin, but de-localized in a 
ring at a radius which depends on the angular momentum 
of the hole. 

In the case of many vortices, simple analytic approxi- 
mations for the states are not known. However, we can 
guess the most important configuration with the help of 
the following arguments [3(|. First we notice that for 
a large number of particles we can replace the center of 
mass in Eq. (|16fl with the origin, i.e. with a fixed point. 
In the same spirit we can assume several localized vor- 
tices at fixed points. When the number of vortices is 
small, their geometrical arrangement will be on a ring. 
In this case the wave function corresponding to Eq. Hlbj) 
would be 



N N 

ae iai ) x ••■ x T\(zi, -ae lak )-$ 



ii 

N 

H(z!;-a n )9 MDD 



IK' 

3k 



'MDD 



(18) 
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where k is the number of vortices, a is distance of the vor- 
tices from the origin and cxj — 2irj/k. Clearly, the above 
wave function does not have a good angular momentum. 
Projecting to a good angular momentum means collect- 
ing out states with a given power of a. We get a state 



*M„ DD +M 



Pm - *i) CX P (- E l Zfe ! 2 /^ (21) 



and for bosons, 



K 



kV 



k(N- 



MDD 



(19) 



which now corresponds to a good angular momentum 
M = Mmdd + kK. Expanding the polynomial in 
Eq. 119(1 shows that the most important configuration 
has the form 



P M exp -El z fclV4 



(22) 



The only difference between the boson and fermion states 
is then the product that makes the fermion state anti- 
symmetric. (In the case of the Laughlin state with q > 3 
the symmetric polynomial is simply P( q -x)N(N-l)/2 — 



no* 



1 1111000 1111111111 0000000) , (20) 

N — K K 

where the number of adjacent zeros between the one's 
equals the number of vortices, and the number of ones 
after the zeros equals to K. 

In the case of bosons, a trial wave function can be 
constructed in a similar fashion, the only difference being 
that the MDD wave function, Eq. 1(15(1 . is replaced by the 
Bose condensate, exp(— \zt\ 2 )- This means that the 
trial wave functions for bosons and fermions are similar 
apart from the product in Eq. 115(1 . which makes the 
fermionic state antisymmetric. 

Though originally developed to describe fractional 
quantum Hall states, a tool that has shown useful to 
describe the high angular momentum states of rotating 
bose condensates |25l l38l l65l Irjtl loll lo8| is the compos- 
ite fermion model of Jain [69j. In this method, the wave 
function for any M is determined first using also higher 
Landau levels and then projected to the lowest Landau 
level. While it would be expected to be applicable in the 
regime M ~ TV 2 , it was more recently extended success- 
fully to describe angular momenta even before the unit 
vortex |65ll70| . 

The construction of the composite fermion state is 
complicated due to the projection to the LLL and does 
not easily reveal the nature of the solution. Neverthe- 
less, the composite fermion picture shows that there is a 
relation between the states of the spinless fermions and 
bosons. In fact, it has been shown [7l| that in the an- 
alytically solvable model of harmonic interparticle inter- 
actions, this relation between the bosonic and fermionic 
solutions is exact and the states can be written as a prod- 
uct of a symmetric homomorphic polynomial Pm and the 
boson or fermion condensate for fermions, 



V. CONCLUSIONS 

At extreme angular momenta, a few particles in a har- 
monic trap crystallize, indpendent of their fermionic or 
bosonic nature. In much analogy to what was observed 
in quantum dots at strong magnetic fields 0, [55[ , this 
crystallization is apparent from the regular oscillations in 
the quantum many-body spectrum: Cusps at the yrast 
line, and their periodicity in angular momentum, are eas- 
ily understood in terms of the simple geometries of these 
so-called Wigner molecules, being of either fermionic or 
bosonic nature [U ■ 

At moderate angular momenta, however, it is known 
that vortices may form, independent of the system be- 
ing bosonic or fermionic. In a harmonic confinement, 
these "holes" that are penetrating the rotating quantum 
system, arrange in simple polygonal structures, very sim- 
ilar to those of the localized Wigner crystallites of parti- 
cles. They lead to similar cusp states in the yrast line, 
with oscillation periods simply following the underlying 
symmetry of the "vortex crystal" . We showed how this 
mapping could be understood in terms of particle-hole 
duality, holding for the bosonic as well as the fermionic 
case |29|. 
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